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Abstract 



We study the 2D vortex-free XY model in a random field, a model for 
randomly pinned flux lines in a plane. We construct controlled RG recursion 
relations which allow for replica symmetry breaking (RSB). The fixed point 
previously found by Cardy and Ostlund in the glass phase T < T c is unstable 
to RSB. The susceptibility x associated to infinitesimal RSB perturbation in 
the high-temperature phase is found to diverge as \ cx (T — T c )~ 7 when T — > 
T 6 + . This provides analytical evidence that RSB occurs in finite dimensional 
models. The physical consequences for the glass phase are discussed. 
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One of the most fruitful methods to study disordered systems is the replica method 
H which allows to average over disorder by introducing n coupled copies of the system. 
Since it reduces the problem to a translationally invariant system this method has the 
considerable advantage, in principle, to allow the use of standard field theory techniques. 
However the limit n —>■ has to be taken, a rather subtle procedure. Indeed, as is well 
known from the large body of works on spin glasses || the free energy functional is in general 
minimized, in a glass phase, by a solution which breaks replica symmetry spontaneously, in 
the hierarchical manner discovered by Parisi ||. Up to now, however, spontaneous replica 
symmetry breaking (RSB) has been clearly demonstrated only when a saddle point method 
could be used naturally. This is the case either in infinite-range or mean field models 
where the saddle point method provides the exact solution of the problem and RSB is 
ubiquitous ||, or, as found more recently 0, in field theories using a Hartree or Gaussian 
variational method (GVM) which becomes exact in the limit of large number of components, 
in some as yet unelucidated sense. As explicitly constructed in both cases the physics of 
RSB corresponds to systems breaking up in many "states" || . Extensions of RSB to finite 
dimensional systems for which such mean-field theory is not exact have been attempted || . 
However, it is far from obvious, as argued in ||, that it is relevant in real, i.e low dimensional 
physical systems. 

Another extensively tried method for these problems is the renormalization group (RG). 
There replicas have been used so far simply as a graph counting trick, to eliminate dis- 
connected graphs in perturbation theory [[Hj. This procedure, however, implicitly assumes 
replica symmetry before taking the limit n — > 0. Although the RG treats fluctuations ex- 
actly and is therefore more accurate than mean field theory, there is a risk that it will miss 
the physics associated to RSB. Thus one would like to link the two methods. 

A good model to look for a renormalization group that allows RSB is the two dimensional 
XY model in a random field and without vortices (i.e with spins S(x) = e 1 ^^ where tp(x) 
is real- valued in [—00,00]). This model is particularly interesting because it is one of the 
simplest and non trivial models of a "glass" in finite dimension, to which several analytical 
methods can be applied ||lTHT6|l- It is relevant to describe several physical disordered systems, 



such as randomly pinned flux arrays confined to a plane OJL8L15], the surface of crystals 



with quenched bulk or substrate disorder |19|| , planar Josephson junctions [P0|| , domain walls 
in incommensurate solids. In a pioneering work ||12|| , Cardy and Ostlund (CO) studied this 
model using the renormalization group. They used the replica method and introduced n 
coupled XY models which they mapped onto a Coulomb gas with n(n — l)/2 types of vector 
charges. They constructed the corresponding RG equations and performed the limit n — > 
in these equations, thus implicitly assuming replica symmetry. The resulting RG recursion 
relations, valid near T c , possess a non-trivial perturbative fixed point for T < T c at weak 
disorder strength g = g* oc T c — T. CO concluded that a "glass" phase exists, controlled by 
this new fixed point. In this phase, one coupling constant is found to runaway to infinity, a 
somewhat peculiar feature. These results were extended in [|l^,^l|], and the disorder averaged 
correlation function C(x) = (ip(x) — ip(0)) 2 was found to grow as C(x) ~ 5 (log \x\) 2 , faster 
than C(x) ~ Tlog \x\ which holds in the high temperature phase and for the pure system. 

Although there is presently agreement that a glass phase exists in this model for T < 
T c its physical properties remain controversial, despite the large number of analytical and 
numerical studies. Two recent numerical simulations [P^,p3] have shown somewhat puzzling 
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results which shed some doubt on the Cardy Ostlund RG. Both simulations were found 
incompatible with the C(x) ~ -B(log \x\) 2 prediction. In the dynamics |22 it was found that 



the velocity develops a non linear dependence on the driving force for T < T c , providing 
evidence of a glassy phase, but with an exponent inconsistent with the prediction of a 
conventional dynamic RG study [ |D| . 

These simulations, however, seem compatible with our previous analytical results using 
the Gaussian variational method [|T^,|T^]. We found that a one-step replica symmetry broken 
solution described the glass phase T < T c with C(x) ~ T c log \ x\ at variance with the result 
of CO, while C(x) ~ Tlog|x| for T > T c . Such a discontinuity in the slope was observed 
in |23[. The simulation |22[ was performed at such weak disorder that the system size was 
shorter than the length predicted by the GVM Jl5[ beyond which glassy behavior can be 
observed in static quantities. One can argue that a variational method, such as the GVM, 
is approximate even in pure systems and misses some of the effects of non linearities better 
captured by the renormalization group. On the other hand, this solution contains the feature 
of RSB and its compatibility with numerical calculations suggests that it has some relevance 
for the physics. Thus, it would be quite interesting to have a direct evidence, in the present 
model, that spontaneous RSB do occur. 

In this paper we present such an evidence. We construct renormalization group recursion 
relations, for the random field XY model, for couplings between replicas of arbitrary sym- 
metry. It contains the Cardy Ostlund recursion relations as a particular case. It allows for 
a more general way of taking the limit n — > 0, using Parisi's type matrices, while remaining 
in the perturbative regime of the renormalization group. The new operators introduced in 
the limit n — > are marginally relevant at T = T c and thus can be incorporated in the 
framework of RG. We show that the low temperature phase does exhibit spontaneous RSB. 
In the high-temperature phase where the RG flows to weak coupling and is therefore cer- 
tainly correct, we compute the linear response to a small RSB perturbation. The associated 
susceptibility diverges at the transition as \ oc (T — T c )~ 7 when T — > T c + . For T < T c the 
replica symmetric flow and the CO fixed point, are unstable to a small RSB. This is to our 
knowledge the first physical model where this effect can be demonstrated in a controlled 
way. 

The Hamiltonian of the 2D vortex-free XY model in a field of random amplitude and 
direction, reads: 



H = J d 2 x |(V^(a;)) 2 - r)(x) ■ Vif){x) (1) 
— Ci(^) cos(^(a;)) — (2(%) sin(if>(x)) 



with r]i(x)rij(x') = Ao(%5(x — x') and d(x)Q(x') = A5ij5(x — x') are two Gaussian white 
noises. This model also describes flux lines with displacements u and average spacing a in 
d = 1 + 1 dimensions, with if) = 2nu/a. A is proportional to the amplitude of disorder with 
Fourier component close to 27i/a and A is the long wavelength disorder. Note that A is 
generated if not present originally in the model. 

After replication of ([!]) and averaging one obtains: 



H. 



off 



J d 2 xj2 

ab 



^V0 a V0 fe - cos(2(<n*) - 



(2) 



3 



where we set g aa = 0. We have used for convenience: 

! 4?rc _ 47rAp g ab _A_ 

and <f) a = ip a /2. This defines the "bare" or starting values parameters, which of course 
are replica symmetric. We now consider a more general situation where the "renormalized" 
parameters K ab and g ab have arbitrary symmetry with respect to the group of replica per- 
mutation. We use the parameterization: 

K ab = 5 ab ~ Kb (4) 

One defines the connected part k c = J2b kab- As will be obvious later there is a transition at 
temperature T c = 4nc in the model ([I]). Using (|3|) and (Q) this corresponds to k c = 0, and 
one has more generally k c = T ~ Tc = r, where r is the reduced temperature. 

In previous applications of RG to disordered systems, recursion relations are established 
for arbitrary number n of replica. Then replica symmetry is assumed and the RG equa- 
tions become simple functions of n. The continuation n — > is then easily taken. In that 
respect the use of replica is a trick of graph counting. In fact one can generally establish 
identical RG equations directly by considering disorder propagators, a method which we 
call "replica symmetric perturbation theory". However, in a glassy system, where many 
metastable states exist, one can question the validity of such a renormalization group proce- 
dure. A crucial assumption of the RG is that one can simply integrate over the short scale 
degrees of freedom independently of the larger scales, to produce an equivalent renormalized 
Hamiltonian. It is not obvious that such a separation between scales exists in glassy systems. 
Short scale degrees of freedom may well be determined by the local minimum they belong 
to, thus depend in effect on larger scales. The presence of long range effects is supported 
by the massless modes found in the expansions around Parisi's solution M. However, the 
replicated Hamiltonian is translationally invariant and all these problems are buried in the 
proper taking of the limit n —>■ 0. In the case where the standard replica symmetric limit 
n — > does not work, a Parisi type RSB might allow to construct a correct RG for glassy 
systems. We will therefore construct RG equations keeping the full matrix structure of the 
couplings. Since we are looking for a scale invariant theory near the transition, we consider 
perturbation theory in r, in the small matrix parameters 2k ab — k aa — k bb and in g ab . This 
leaves unconstrained the degree of freedom of a constant shift of all the elements of k ab . K~ b 
is easily inverted as K ab = (1 — k c )~ 2 [{\ — 2k c )S ab + k ab ] ~ 5 ab + k ab , which is valid even if 
the k ab are large provided the above-mentioned parameters are small. 

The constants g ab , k ab and k c = r correspond to quadratic interactions in an equivalent 
fermion problem |24| . Using standard renormalization group methods, either on the fermion 



form or on the classical Hamiltonian (§) p4j , one can derive to second order the general RG 
for the replicated system 

-^j- = (2& a 6 — k aa — k bb )g ab + - g ac gbc 
dk ab 1 2 dr 

= i 9 *> n = (5) 
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where I is the standard logarithmic scale. In these equations a ^ b is implied. Note that the 
reduced temperature r is unrenormalized. The subspace of replica-symmetric parameters 
9ab = 9, k a ^b = k is obviously preserved by the RG and the equations (|5|) for g, k and r 
reduce to the one obtained by Cardy and Ostlund in the limit n = 0: 

In the high temperature phase, r > 0, the disorder g is irrelevant and the fixed point is a pure 
Gaussian system with C(x) = ((f>(x) — </>(0)) 2 ~ (|r + 0(g)) log \x\. The renormalization of 
the off-diagonal part k contributes to 0(g) since k goes to a constant (which stays finite up to 
t = 0). In the low temperature phase, there is a fixed point at g* = 2|r|. Note however that 
flows to infinity. This is a peculiar situation which within the replica symmetric scheme 
does not lead to inconsistencies, since k does not feedback to any order in perturbation 
theory (only averages (Z)aO a a ) 2 with J2 a C a = appear). The flow of k was used [[HJ to 
predict C(x) ~ £>(log |x|) 2 . 

Let us now look at the RG flow for coupling constants parameterized by Parisi matrices. 
g a b and k ab , a ^ b are replaced by functions g(u), k(u), with < u < 1. The symmetric case 
corresponds to constant functions. The independent variables are now r, g(u) and k(u). 
The equations now read: 

d9{U) -' -2r + 2(k(u) - (k)))g(u) - (g)g(u) - ~ f dt^i,) - g{u))\ = -a{uf 



dl K ' ' 2 Jo v v y v " dZ 4' 

(7) 

where (fc) = Jq 1 dvkiy). These functions depend on I but this will be written explicitly only 
when needed. Some features are immediately apparent on (0). There is now a spectrum of 
dimensions, given by — r + k{u) — (k), for the operators corresponding to the non replica 
symmetric couplings g a b- The dimension of these operators should be small to neglect the 
effect of higher order replica terms. Also k(u) now feeds back in the equation for g(u), thus 
one anticipates an instability. 

Let us first check the stability of the replica symmetric flow in the low temperature phase. 
We separate the replica symmetric part and write g(u) = gi + e(u) and k{u) = ki+A(u) with 
(e) = and (A) = 0. Starting from the CO fixed point g = g*, ki being the running coupling 
constant when g = g*, one obtains from (|7|) the deviations from the replica symmetric 
solution to linear order: 

deju) n * A , s dAju) g* 

Thus the RG flow is clearly unstable to replica symmetry breaking when r < 0. The 
eigenvalue of instability is A = 2|r|. When r < there is no small coupling fixed point and 
the flow goes to strong coupling [25]. This is clear from (|7|) since A(u) can only reach a 



fixed point if g 2 (u) — (g 2 ) goes to zero. This is a strong indication that the low temperature 
phase corresponds to a replica symmetry broken solution. To conclude unambiguously on 
that issue, we now consider the high temperature phase r > where the RG is exact since 
there is a weak coupling fixed point g(u) = (and A(u) small). By computing the linear 
susceptibility to small RSB perturbation we show that spontaneous RSB occurs at T c . 
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Let us define the susceptibility to RSB as follows. We start, for r > with a given 
disorder g and add a small RSB perturbation g(u) = go + eo( M ) uniform in space. This 
will result in replica non symmetric part in the correlation functions (0 a (<?)0&( — <?)) = {^ab + 
k* ab )/q 2 where k* b are the renormalized coupling at the fixed point g{u) = 0. A susceptibility 
is defined as: 



A*{u) = k*(u) - (k*) = X e (u) (9) 

in the limit eo — > 0. One should consider a \ function of u but it turns out that it has the 
above simple form. Alternatively, one can add a quadratic RSB perturbation 5fc a t,V0 a V0fo 
and define a corresponding susceptibility yf as in (||). To compute \ let us expand the RG 
equations (|7]) to linear order in e . As can be seen from (0), A can be considered as linear 
in e . 

dg 2 dA(u) 1 , , 

Ti = ~ 2t9 ~ 9 > — = 2 9e{u) 

^ = -(2 T + g)e(u) + 2gA(u) (10) 

At linear order, each u can be treated independently. We perform the following rescalings 
x = 2rl, e{u) = e (u)f, A{u) = €o(u)a, g = 2rg. The first equation of (|T0|) integrates to give 
g( x ) = a+ i-e-* wriere we have defined a = 2r/g . Defining Z(x) = | J X dyA(y) + 1, ([TO]) 
gives f(x) = ae~ x Z(x)/(a + 1 — e~ x ) where Z(x) satisfies: 

= Z(x)e~ 2x /(a + 1 - e~ x ) 2 (11) 

with initial conditions, Z(0) = 1 and Z'(0) = 0. The problem depends thus only on one 
variable a = 2r/g . The susceptibility is obtained from the asymptotic value A^ and is 
X = \im x _ QO aZ'(x)/2. 

A full solution of (|TT| ) can be obtained numerically, but \ can be estimated from con- 
sidering the two asymptotic regimes a ^> 1 and a <C 1. In the first regime a ^> 1, 
(|TT|) becomes Z"(x) = e~ 2x Z(x)/a 2 , whose solution is in terms of Bessel functions, 
Z[x) = Cil (e~ x / a) + c 2 K (e~ x / a) . Ci j2 are two constants determined from the initial 
conditions. This gives in that regime: 

X=-h(-)~- (12) 

To investigate the behavior close to the transition one has to look at a 1 and to match 
two regimes in x. For i < 1, (11) can be written Z"[x) = Z[x) /{a + x) 2 whose solution are 



power laws Z oc (a + x) v where v satisfies the golden mean equation v 2 + v = 1. We define 
7 = (y/E - i)/2 ~ 0.618. The solution is Z(x) = ((1 + x/a)" 7 + 7 2 (1 + x/a) 1 / 7 )/(l + 7 2 ). 
When x becomes of order 1, one can use the equation of the first regime with a — *■ a + 1. A 
simple matching at x ~ 1 gives: 

X oc C(g /2ry (13) 
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We performed a numerical integration of (|TT|) which confirms the analytic estimates (|T2l , |T3|) 
and gives C ~ 0.165. x' can be computed similarly from pi]) , with modified initial conditions 
Z(0) = and Z'(0) = 2/ a and has the same divergence as x- There is also a non linear 
response in eo in the high temperature phase. A runaway flow occurs for smaller and smaller 
values of eo when r — ► 0, roughly when x e o > T - 

The divergence of x m dUP when r — ► fixed go shows that RSB occurs spontaneously 
|25fl for t < 0. Quantities like (4>a(<l) < fib(—q)) acquire a replica non symmetric part. Defining 
Qab = (V0 a (x)V0 fe (x)), a possible order parameter for RSB is Q ab = 2q ab +J2 c ^ a Qac+Hc^b Ibc- 
The perturbation Sk a b is the conjugated field to Q a b- Thus x' is the inverse effective mass of 
Q a b, and its divergence signals an instability a la de Almeida Thouless |?7| when T < T c . By 
analogy with the sine Gordon model, a reasonable interpretation of the runaway flow of g a b 
is that a "mass" develops for some modes below T c . This scenario is in good agreement with 
the findings of the variational method [|n|[l||, where correlation functions become non replica 
symmetric and a fraction 1 — u c of the modes become massive below T c , with u c ~ T/T c . 
These modes correspond, in the Coulomb gas, to types of charges which unbind. Note that 
a fraction u c remains massless (corresponding to bound types of charges). This method, 
being non perturbative contrary to RG, allows for the generation of a mass. 

In this paper, we have shown that replica symmetry breaking can be incorporated in the 
RG. We found new relevant operators which break replica symmetry when n — ► for T < T c . 
The associated susceptibility diverges in the high temperature phase T > T c . The replica 
symmetric fixed point found by Cardy Ostlund is unstable. Observable consequences of RSB 
could be looked for in the dynamics of this system. RSB is usually accompanied by effects 
of aging, persistent correlations, and breaking of fluctuation dissipation theorem p8| . In 
our opinion it is also probable that the static correlation functions should be different from 
the one predicted by the symmetric RG, as hinted at by the variational method. Further 
numerical and experimental results would be of great interest. 

We are grateful to E. Brezin, M. Gabay, T. Hwa and Y. Shapir for useful discussions. 
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